ASYMPTOTIC STABILITY OF PERIODIC SOLUTIONS FOR 
NONSMOOTH DIFFERENTIAL EQUATIONS WITH APPLICATION 
TO THE NONSMOOTH VAN DER POL OSCILLATOR 
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Abstract. In this paper we study the existence, uniqueness and asymptotic stability of the 
periodic solutions for the Lipschitz system x = £g{t, x, e). Classical hypotheses in the periodic case of 
second Bogolyubov's theorem imply our ones. By means of the results established we construct, for 
small e, the curves of dependence of the amplitude of asymptotically stable 27r— periodic solutions of 
the nonsmooth van der Pol oscillator il+e {\u\ — 1) u+{l + ae)u = eAsint, on the detuning parameter 
a and the amplitude of the perturbation A. After, we compare the resonance curves obtained, with 
the resonance curves of the classical van der Pol oscillator ii + e — l) it + (1 + ae)u = eAsint, 
which were first constructed by Andronov and Witt. 
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1. Introduction. In the present paper we study the existence, uniqueness and 
asymptotic stabihty of the T-periodic solutions for the system 

(1.1) x = eg{t,x,e), 

where £ > is a small parameter and the function g e C'^(M x R*^ x [0, 1],]R''') is 
T-periodic in the first variable and locally Lipschitz with respect to the second one. 
As usual a key role will be played by the averaging function 



(1.2) goiv) = J giT,v,0)dT, 



and we shall look for those periodic solutions that starts near some vq G gg"^(0). 

In the case that g is of class C^, we remind the periodic case of the second 
Bogolyubov's theorem ([6], Ch. 1, § 5, Theorem II) which represents a part of the 
averaging principle: det (go)'(i^o) 7^ assures the existence and uniqueness, for e > 
small, of a T-periodic solution of system lll.l]) in a neighborhood of vq, while the 
fact that all the eigenvalues of the Jacobian matrix (go) '(I'd) have negative real part, 
provides also its asymptotic stability. This theorem has a long history and it includes 
results by Fatou , Mandelstam-Papaleksi |3D] and Krylov-Bogolyubov 251 § 2] . 

Second Bogolubov's theorem gave a theoretical justification of resonance phe- 
nomenons in many real physical systems. The most significant example is the classical 
lamp oscillator whose scheme is drawn at Fig. 11.11 and whose current u is described 
by the following second order differential equation 

(1.3) il+-^{RC- Mi'{u)) ii + u}^u = -^F{t), 
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Fig. 1.1. Circuit scheme for the classical lamp oscillator (see Andronov- Vitt- Khaikin fllj. 
Ch.VIII, §a, Fig. 348, Malkin fSSS, Ch.I, §5, Fig. 1, Nayfeh-Mook fSSS, ^3.1.7, Fig. 3-5). 




Fig. 1.2. Characteristics of the lamp of the circuit of Fig. (a) - lamp in a harsh regime 

(see Andronov- Vitt- Khaikin Ch.IV, §7, Fig. 212b, Malkin ^gSj, Ch.I, §5, comments for Eq. 5.3- 
5.4); (b) - lamp with saturation (see Andronov- Vitt- Khaikin fTV . Ch.VIII, §5, Fig. 364)', (<^) - lamp 
without saturation (see Andronov- Vitt- Khaikin fT]', Ch.IX, §7, Fig. 482) 



of periodic solutions in this system is performed in almost every book on nonlin- 
ear oscillations (see Andronov- Vitt-Khaikin [T], Ch.VIII, §2, Malkin [2S], Ch.I, §5, 
Nayfeh-Mook [35], §3.1.7) but with approximation i{u) — i(a){u) = So + Siu — ^S^u^ 
(leading to the classical van der Pol equation). Therefore, it is natural to look for 
a technique that permits avoiding this approximation and dealing with the original 
shape of the lamp characteristic drawn at Fig |1.2b . that expects to give more accurate 
correspondence between theoretical and experimental results. Moreover, a wide class 
of physical systems is modelled by circuit Fig ll. II whose lamp either has or has not a 
saturation that leads to the characteristic drawn at Figures [TT^ and ll.2b respectively. 
Though the unforced equation (|1.3p (i.e. for F = 0) with i described by Fig. 11.2b 
and Fig.[Ll; is weU studied (see Andronov- Vitt-Khaikin [J, Ch.VIII, §3 and Ch.IX, 
§7 respectively), the question about resonances in these equations when F ^ (e.g. 
F{t) = eXsmt) is still open. With regard to equation (|1.3[) with lamp characteristic 
given by Fig. ll.2b . ll.2b or ll.2b we finally note that Levinson's change of variables (see 
|28j . pass from Eq. 2.0 to Eq. 2.1) allows to rewrite equation (|1.3p as the following 
system 

ziRC + M{i{zi)~i{0)) ziRoC + Mo{i{zi) ^ i{0)) 

Zl — Z2 — Z2 — e- 



LC LC 
1 



Z2 = —ujzi -\- F{t) = —zi -\- £ ( hzi -\ — Asint 



whose solution (zi, Z2) gives a solution u = zi to (|1.3p . Then the change of variables 
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brings this system to the form (|l.ip with 



\ _ f cost sint \ / xi{t) 
Z2{t) ) \ — sini cost j \ X2{t) 
Lipschitzian (in the second variable) g. Therefore the goal of the paper is to generalize 
second Bogolubov's theorem for the case when g in p.ip is Lipschitzian. 



L 




Fxii) 



Fig. 1.3. Driven Chua's circuit (see 03', W3, \M' 

Another motivation of this paper comes from modern electrical engineering where 
driven (or forced) Chua's circuit drawn at Fig. II. 31 is a subject of an enormous number 
of papers. Circuit at Fig. 11.31 is described by the following three-dimensional system 



'-'i—n- = H^'Ci) + Fi{t), 



(1.4) 



dt R 

dvc2 _ - 
dt ^ 
diL 
~dt 



R 



L 



-VC2 + F2{t,vc2) 



where i{v) is the characteristic of the Chua's diode whose shape drawn at Fig. 11.41 is 
piecewise linear. Many numerical simulations have been suggested around dynamics 




Fig. 1.4. Nonlinear characteristic of the Chua's diode of the circuit drawn at Fig. U.31 given 
by i{v) = GbV + (l/2)(Ga — Gh) {\v + Bp\ — \v — Bp\) , where Ga, Gf,, Bp S R are some constants 
depending on the properties of the Chua's diode (see Chua 



of ([SI) in the recent literature, see [39], [20] for Fi 7^ and Fa 7^ 0, [5], |34] for 
Fi = and periodic F2, [11] where both Fi and F2 are periodic. Generalization 
of the second Bogolubov's theorem for equations (|l.ip with Lipschitzian right hand 
part will allow for the first time theoretical detection of asymptotically stable periodic 
solutions in certain configurations of the driven Chua's circuit (|1.4p provided that Ci 
is large enough. This promises to forestall some numerical simulations (e.g. to work 
out interesting parameters of the driven Chua's circuit) giving a significant impact 
for further experiments. 

For a large extent the phenomenal interest in generalizing of the second Bogol- 
ubov's theorem for the Lipschitzian case comes from mechanics, where systems with 
piecewise linear stiffness describe various oscillating processes. A prototypic process 
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Fig. 1.5. A prototypic device (a) where a driven mass is attached to a immovable beam via a 
spring with piecewise linear stiffness (b), see e.g. JSjj, \24i (Ch.I, p. 16 and Ch.IV, p. 100) and \37^ . 



of this type is exhibited by the device drawn at Fig. 11.5b . where a forced mass is at- 
tached to a spring whose stiffness changes from ki to ki+k2 when the mass coordinate 
crosses in the negative direction. This device is governed by the folfowing second 
order differential equation 

(1.5) mx + P{x) ^ F{t,x,x), 

where piecewise hnear stiffness P is drawn at Fig. 11.5b . Depending on a particu- 
lar configuration of the device of Fig. 11.5b various terms can stay for F in (|1.5[) . 
It is F{t,x,x) = —f{x)x + McosLot with piecewise constant / for shock-absorber 
and jigging conveyor (see [23], Ch.I, p. 16 and Ch.IV, p. 100 where original second 
Bogolubov's theorem is employed without justification). Levinson's change of vari- 
ables (see transformation of Eq. p.3p above) allows to rewrite (|1.5p as a Lipschitz 
system. It takes simpler form F{t, x, x) = —cx + MQ{t) for an impact resonator and 
F{t, X, x) — —cx + M sin for a cracked-body model (see |37j and [8], where numer- 
ical experiments are performed solely). In each of these situations equation (|1.5p can 
be rewritten as system (jl.ll) with Lipschitzian g provided that the constant ^2 and 
the amplitude of the force F are sufficiently small. Therefore the related generaliza- 
tion of the second Bogolubov's theorem promises to justify or explain the resonances 
appeared in the cited results. We note that the recent report by Los Alamos Na- 
tional Laboratory [14] describes increasing interest in a specific form of the model of 
Fig. 11.5b . called cracked-body model and, particularly, in suspension bridges models 
that is why we discuss the contribution of the present paper to the latter model in 
details. 

The first idealization of a one-dimensional suspended bridge is drawn at Fig. 11.6b . 
It is represented (see [T7], [IS]) by the beam bending under its own weight and sup- 
ported by cables whose restoring force due to elasticity is proportional to u"*" (see 
Fig. 11.6b ). where u — u{t,x) is the displacement at a point at distance x from one 
end of the bridge at time t and u measured in the downward direction. Looking for 
u of the form u{t,x) = z{t) sm{Trx/ L) and considering F[x,t) = h{t) smijrx / L) we 
arrive (see [T^) to the following form of equation (|1.5p 

(1.6) mz + 6z + c{7r/L)^z + dz^ ~ mg + h(t), 

where the constant m > is the mass per unit length, 5 > is a small viscous 
damping coefficient, c > measured the flexibility or stiffness of the bridge, L > 
is the length of the bridge, d > represents the stiffness of nonlinear springs and h 
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Fig. 1.6. (a) - The first idealization of the suspension bridge: the beam bending under its own 
weight is supported by the nonlinear cables (see Lazer-McKenna }26f . Fig. 2); (b) - characteristic 
of stiffness of nonlinear springs. 



is a continuous T-periodic force modelling wind, marching troops or cattle (see |19| 
for details). Considering c > and d > fixed and assuming that either c > and 
h{t) are sufficiently small, or c > is fixed and h{t) is sufficiently large, or c > is 
sufficiently small and h{t) fixed, Giover, Lazer, McKenna, Fabry (see [17], [27], [26] . 
[TS] ) proved various theorems on location of asymptotically stable T-periodic solutions 
in (|1.6|) . The question what happens with these solutions when d > appears to be 
the same magnitude small as 5 > and /i(t), is open for a while and can be resolved 
by means of the generalization of the second Bogolubov's theorem we propose. Lazer 
and McKenna proved in [57] that the Poincare map for (|1.6p is differentiable, but 
we note that it is not an enough argument to apply the original second Bogolubov's 
theorem since it requires the differentiability of functions participating in (|1.6p as well. 

At the end of the applications review we note that system (|1.4p describing the 
Chua's circuit (Fig. II. 3p appeared in the recent time to govern mechanical systems 
with so-called "negative slope" (see Awrejcewicz 0], §8.2.2). Thus, analogous to 
Chua's circuit applications of the result of this paper are also possible to these me- 
chanical systems. 

It was Mitropol'skii who first noticed that various applications require generaliza- 
tion of the second Bogolubov's theorem for Lipschitz right hand parts. Assuming that 
g is Lipschitz, 50 G C^{M.^ , M'"') and that all the eigenvalues of the matrix (go)'(wo) have 
negative real part Mitropol'skii developed the second Bogolyubov's Theorem proving 
the existence and uniqueness of a T-periodic solution of system (|l.ip in a neighbor- 
hood of Vq. There was a great progress weakening the assumptions of Mitropol'skii 
in his existence result (see Samoylenko |38j and Mawhin [31 J, but not of his unique- 
ness result. Moreover, the asymptotic stability conclusion of the second Bogolyubov's 
Theorem remained to be not generalized for Mitropol'skii's settings (namely, when g 
is Lipschitz) for a long time. It has been done recently by Buica-Daniilidis in f9l for 
a class of functions v ^ g{t,v,0) differentiable at vq for almost any t £ [0,T], but it 
is assumed in [9] that the eigenvectors of the matrix (5o)'(wo) are orthogonal. 

In the next section of the paper assuming that g is piecewise differentiable in the 
second variable we show in Theorem l2.5l that Mitropol'skii's conditions imply not only 
uniqueness, but also asymptotic stability of a T-periodic solution of system (jl.ip in 
a neighborhood of vq. In other words we show that Bogolyubov's theorem formulated 
above is valid when g is not necessary . Theorem 12.51 follows from our even more 
general Theorem 12 . 1 1 whose hypotheses do not use any differentiability neither of g nor 
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of go- In Section 3 we illustrate our result constructing resonance curves of nonsmooth 
van der Pol oscillator [TB] . This application has been chosen since it allows to compare 
the issues of our Theorem 12.51 with the classical results [5] and [5] by Andronov and 
Witt obtained for original van der Pol oscillator. 

2. Main results. Throughout the paper fl C M*"' is some open set. For any 
5 > we denote Bs{vq) — {v £ R'' : \\v — vo\\ < . We have the following main 
result. 

Theorem 2.1. Let g eC"iRxnx [0,1], R'=) and vq e fl. A ssume the following 
four conditions. 

(i) For some L > we have that \\g(t, vi, e) — g{t, V2, s)\\ < L \\vi — V2\\ for any 

t e [0,r], VI, V2 en, ee [0,i]. 

(a) For any 7 > there exists 6 > such that 

Jo 5(t, f 1 + u{t), e)dT - /J" .g(r, V2 + u(t), e)dT 

- Jo 9iT,vi,0)dT + g{T,V2,0)dT < j\\vi -U2II 

for any u G C°([0, r],R'=), < S, vi,V2 € Bs{vo) and e e [0,(5]. 

(Hi) Let go be the averaging function given by IL2\} and consider that goivo) — 0. 

(iv) There exist g e [0, 1), a, (5o > and a norm \\-\\o on M.^ such that ||wi + ago{vi) 
-V2 - agoiv2)\\o < q\\vi - V2\\o for any vi,V2 € Bsgivo). 
Then there exists Si > such that for every e G (0, ^i] system il.ll) has exactly one T- 
periodic solution x^ with x^{Q) e Bs-^{vo). Moreover the solution x^ is asymptotically 
stable and Xi,{Q) — *■ Wq as e — > 0. 

When solution x{-,v, e) of system (|l.ip with initial condition a;(0, v,e) — v is well 
defined on [0, T] for any v e Bs^,{vo), the map v ^ x{T,v,e) is well defined and it is 
said to be the Poincare map of system (jl.ip . The proof of existence, uniqueness and 
stability of the T-periodic solutions of system (|l.ip in Theorem 12.11 reduces to the 
study of corresponding properties of the fixed points of this map. 

In order to prove Theorem l2.1l we observe from ()1.1|) that x{T, v, e) can be repre- 
sented as 

T 

x{T,v,e) = V + egsiv), where g^iv) J g{T,x{T,v,e),e)dT, 



and we use the following result which claims that properties (P and ([u]) are also 
applied to in a suitable sense. 

Lemma 2.2. Let g e C°(Rxl]x [0, 1],R'') andSa > be such that Bs^ivo) C n. If 
(OP is satisfied then there exist S G [0, Sq] and Li > such that the map (v, s) ^— > ge{v) 
is well defined and continuous on Bs„{'>^a) x [0,(5] and 

\\ge{vi) - geMW < Li\\vi - V2\\ for any s e [0,S], vi,V2 e Bs„{va). 

If both |ip and are satisfied then for any 7 > there exists 6 € [0, ^o] such that 

\\9e{vi) - ffo(wi) - ge{v2) + go{v2)\\ < l\\vi - V2\\ 

for any vi,V2 G Bs{vo) and e G [0, S]. 

Proof. Using the continuity of the solution of a differential system with respect 
to the initial data and the parameter (see [3S], Ch. 4, § 23, statements G and D), 
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we obtain the existence of Eq > such that x{t,v,e) G il for any t £ [Q,T], v G 

Bso (vo) e E [0, Eq]- Using the Gronwall-Bellman Lemma [13, Ch. II, § 11] from the 

t 

representation x(t,v,e) ^ v + s J g{T,x{T,v,e),e)dT and the property Q we obtain 



\\x{t,vi,e) - x{t,V2,e)\\ < e^^'^\\vi - W2II for aU t e [0,T], vi, V2 £ Bso{vo) and 
e e [0,eo]. Therefore y{t,v,e) = g{T,x{T,v,e),e)dT satisfies the fohowing property 



(2.1) \\yit,vi,e)-y{t,V2,s)\\ <ii||wi-V2|| 

for alH e [0,T], vi, V2 G BsJvq), s G [0, £0] and Li = LTe^°^'^. Since ge{v) 
y{T,v,s) the first part of the lemma has been proven. 

Taking into account that x{t, v,s) = v + sy{t, v, e) we have 

(2.2) y{T, VI, e)- y{T, vi, 0) - y{T, V2, e) + y{T, V2, 0) = Ii{vi,V2,e) + l2{vi,V2, e) 
where 

Ii{vi,V2,e) ^ / [g{T,V2 + ey{T,vi,s),e) - g{T,V2 + ey{T,V2,e),e)]dT 



l2{vi,V2,e)^ [ [{g{T,vi + ey{T,vi,e),e) - g{T,V2 + ey{T,vi,e),e))]dT 

T 

{g{T, vi,0) - g{T, V2,0))dT. 

Since {t,v,e) 1-^ y{t,v,e) is bounded on [0,T] x Bsg{vo) x [0, £0], we have that 
ey{t,v,e) ^ as £ ^ uniformly with respect to t e [0,7] and v G Bsf,{vo). 
Decreasing £0 > 0, if necessary, we get that V2 + ey{t,vi,e) G for any t G [0,T], 
vi,V2 G Bsg{vo), £ G [0,£o]. By assumption jij and relation (|2.ip we obtain that 
||/i(wi,W2,£)|| < T ■ eLLi\\vi - W2II for all e G [0, £0], vi,V2 G Bso{vo). 

We fix 7 > and take S > given by ([n]). Without loss of generality we can 
consider that 6 < minj^o, £0, l/i'^TLLi)}. Therefore assumption juj implies that 
\\l2{vi,V2,e)\\ < (7/2)||vi - t;2|| for any £ G [0,6], vi,V2 G Bs{vo). Substituting 
the obtained estimations for Ii and I2 into (j2.2p we have \\y{T,vi,e) — y{T, vi,0) — 
yiT,V2,e) +yiT, V2,0)\\ < (eTLLi + ^/2)\\vi ^ V2\\ < tII^^i " ^^2!! for any £ e [0,6], 
vi,V2 G Bs{vo). Hence the proof is complete. □ 

Lemma 2.3. Let go : il ^ M.^ satisfy assumption Iwj) with some q G (0, 1), a, 60 > 
and a norm ||-||o onM*-'. Then [\vi + ego{vi) — V2 — £go('y2)|lo !i (1 ~ ^(1 ~ l)/'^) ll^i~ 
V2[\o for any vi,V2 G Bs„{vo) and any e G [0,a]. 

Proof. Indeed, the representation v + sgo{v) = (1 — e/a)v + e/a{v + ago{v)) 
implies that the Lipschitz constant of the function / + ego with respect to the norm 
II • llo is {1 — e/a) + e/aq = 1 — e{l — q)/a. □ 

Proof of Theorem \2.1\ By Lemma [2.21 we have that there exists 61 G [0, (5o] such 

that 

(2.3) ||.ge(fl)-ffo(fl)-<?e("2)+5o(f2)||o < ((1 - (7)/(2a)) ||t;i - 1-2 || 

for any e G [0, i5i], vi,V2 G Bs^[vq). First we prove that there exists £i G [0,6i] such 
that for every e G [0, £i] there exists G Bs^ (vq) such that x{-, Vg, e) is a T~periodic 
solution of (jl.ip by showing that there exists such that x{T, u^, £) — v^. Using (pli|) 
and (pv| we have 

[[v + agoiv) - vq[\o < q[\v - t;o||o for any v G Bg^ (uo)- 
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Therefore we have that the map / + agQ maps Bs-^{vq) into itself. From Lemma 12.21 
we have that there exists Eq > such that the map {v,e) i-^ ffe(i') is well defined 
and continuous on Bsj^{vo) x [0,eo]- We deduce that there exists £i > sufficiently 
small such that, for every e € [0, ei], the map / + ag^ maps Bg^ (vq) into itself as well. 
Therefore, by the Brouwer Theorem (see, for example, [23l Theorem 3.1]) we have 
that Bs-^{vo) contains at least one fixed point of the map / + ag^ for any e € [0,£i]. 
Denote this fixed point by v^- Then we have gdv^) — and x{T,Ve,s) — for any 
e G [0,ei]. 

Now we prove that is the only T-periodic solution of (jl.ip originating 

near a-nd that, moreover, it is asymptotically stable. Knowing that x{T,v,e) = 
V + £ge{v) we write the following identity 

(2.4) x(r, v,e)=v + ego{v) + e ig,{v) - go{v)) . 

Using Lemma 12.31 we have from (|2.3|) and (|2.4p that 

\\xiT,v,,e) - x{T, V2,s)\\^ < (1 - £(1 - q)/a + e(l - q)/{2a))\\vi - V2\\o 

= (l-£(l-g)/(2a))||i;i-z;2||o, 

for all vi,V2 e Bs^{vo) and e G [0, We proved before that there exists ei > 
that, for every e £ [0,ei] there exists G Bs-^{vo) such that a;(-,We,e) is a T-periodic 
solution of (|l.ip . Since e{l — q)/{2a) > and Si < Si the last inequality implies that 
for each e G [0, the T-periodic solution a;(-,We,e) is the only T-periodic solution 
of (jl.ip in Bsi{vo) and, moreover (see [23i Lemma 9.2]) it is asymptotically stable. □ 

Remark 2.4. We note that a similar result close to Theorem \2.1\ is obtained by 
Buica and Daniilidis (see Theorem 3.5). But instead oj the assumption |wi|) with 
fixed a > it is assumed to be satisfied for any a > sufficiently small. Although, 
Lemma \2.S\ now implies that it is the same to assume for only one a > and, 
respectively, for all a > sufficiently small. The advantage of our Theorem \2.1\ is 
that it does not require differentiability of g{t, •,£) at any point, while needs it at 
Vq. See also Remark \2.8\ 

In general it is not easy to check assumptions and ((ivl) in the applications 
of Theorem 12.11 Thus we give also the following theorem based on Theorem 12.11 
which assumes certain type of piecewise differentiability instead of (|Ii]) and deals with 
properties of the matrix {g^)' [vq) instead of the Lipschitz constant of go- 

For any set M C [0, T] measurable in the sense of Lebesgue we denote by mes(Af ) 
the Lebesgue measure of M (see ^ , Ch. V, § 3). 

Theorem 2.5. Let g G C"(M x 17 x [0, 1],M'') satisfy 0). Let go be the averaging 
function given by il.2\) and consider wq G 51 such that go(^^o) = 0. Assume that 

(v) given any 7 > there exist 5 > and M C [0, T] measurable in the sense 
of Lebesgue with mes(M) < 7 such that for every v G i?j(uo), t G [0,T] \ M 
and e G [0,(5] we have that g{t,-,e) is differentiable at v and \\g'^{t,v,e) — 

g'^{t,vo,m<l- 
Finally assume that 

(vi) go is continuously differentiable in a neighborhood of vo and the real parts of 
all the eigenvalues of {go)' {vo) are negative. 

Then there exists > such that for every e G (0, Si], system 11. 1\) has exactly one T- 
periodic solution x^ with Xe{0) G Bs-^{vo). Moreover the solution x^ is asymptotically 
stable and a;e(0) — s- uq as £ — s- 0. 
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For proving Theorem 12.51 we need two preliminary lemmas. 

Lemma 2.6. Let g e C°(M x x [0, IJjR'^) satisfying If holds then ^ is 
satisfied. ^ 

Proof. Let 7 > be an arbitrary number. We show that ^ holds with 5 — 5/2, 
where 6 is given by (jvj) applied with 7 — min{7/(4i), 7/(4r)}. We consider also 
M C [0, T] given by (jvj) applied with the same value of 7. 

Letue C''([0,r],R'=), ||w|| < ^andi^(w) = giT,v + u{T),e)dT~ g{T,v,0)dT. 
Let vi,V2 e Bs{vo) and e G [0,6]. We have F{v) = Fi{v) + F2{v), where Fi{v) = 
iMigiT^v + uir), e)~g{T, v, 0))dT and F2{v) = /j^, ,j,j^j,^(5(t, w + m(t), e)~g{T, v, 0))dr. 
By we have that \\Fi{vi) - Fi{v2)\\ < 2L ■ mes(M)||t;i - W2II < 2L^\\vi - wall < 
(7/2)||t;i — W2II. On the other hand, using (jvj), we will prove that a similar relation holds 
for F2. In order to do this, we denote /i(r, v) — g{T, v + u{t), e) — g{T, v, 0). Notice 
that for each r S [0,T]\M we can write h'^{T, v) — {g'^{T, v + m(t), e) — g(,(T, vq, 0)) — 
{g!u{T, V, 0) — g'^{T, vq, 0)). As a direct consequence of (jvj) we deduce that w)|| < 

27 for all V G Bs{vq) and r G [0, T] \ AI . Now applying the mean value theorem for the 
function /i(t, •), we have ||/i(r,ui)-/i(T,U2)|| < 27||?;i-U2|| for aU r e [0, r]\M and all 
vi,V2e Bs{vq). Then ||F2(?;i)-i^2(v2)|| < / ||/i(t, wi) - /i(t, ?;2)||dr < 2T7||z;i - 

[0,T]\M 

^^211 < (7/2)||ui — W2||- Therefore, we have proved that — ^"("2)11 < 7||wi — W2II, 

that coincides with (jul. Q 

Lemma 2.7. Let g^ : Vl ^ R**' satisfying assumption Iw)) for some vq G ft. Then 
there exist q G [0, 1), a, 60 > and a norm \\ ■ \\q on M'^' such that is satisfied. 

Proof. If A is an eigenvalue of ctigoYivo) then A+1 is an eigenvalue of I+{agoy (vq). 
Since the eigenvalues of a{gQy{vo) tends to as a —> and have negative real parts 
then there exists a G [0, 1) such that the absolute values of all the eigenvalues of 
/ + a{goy{vQ) are less than one. Therefore (see p. 90, Lemma 2.2]) there exist 
g G [0, 1) and a norm || • ||o on R'^ such that sup||^|||^<i \\£_ + ck(go)'(wo)^||o < 9- 

By continuous differentiability of ga in a neighborhood of vq we have that \\ga{vi) — 
9o{v2) - (.9o)'(wo)(wi - V2)\\ /\\vi - W2II < \\9o{vi) - 30(^2) - {goYMivi - W2)|| + 
WigoY {v2){vi-V2)- (goY M{vi-V2)\\ /\\vi-V2\\ ^ as max{ | j wi - || , || W2 " ^0 1| } ^ 0. 
Therefore taking into account that all norms on R'^ are equivalent, there exists > 
such that \\go{vi) - .90(^2) - {goYM{vi - V2)\\o < (1 - 9)/(2a) - W2II0 for all 
vi,V2 e Bsoivo). Then 

\\vi+ago{vi) -V2- ago{v2)\\o 

< a\\go{vi) - go{v2) - {goYMivi - W2)||o + - V2 + a{goY{va){vi - V2)\\o 

< {l + q)/2\\vi-V2\\o, 

for all vi,V2 G Bso{vo). □ 

Proof of Theorem \2.5[ Lemmas 12.61 and 12.71 imply that assumptions ^ and (pv|) 
of Theorem l2.1l are satisfied. Therefore the conclusion of the theorem follows applying 
Theorem Em □ 

It was observed by Mitropol'skii in [33] that in spite of the fact that g{t, •,£) in 
(jl.ip is only Lipschitz, function go turns out to be differentiable in applications. In 
particular, one will see in Section 3 that this is the case for the nonsmooth van der 
Pol oscillator. 

Clearly if g G C^R x R'' x [0, 1],R'=) then (P and (jv} hold in any open bounded 
set ft C R'^'. Therefore Theorem 12.51 is a generalization of the periodic case of the 
second Bogolyubov's theorem formulated in the introduction. 

Remark 2.8. Our Theorem \2.5\ does not require that the eigenvectors of{gQY{vo) 
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be orthogonal as in the result of Buica and Daniilidis ([9^, Theorem 3.6). Moreover 
assumption (H2) of 19] is more restrictive than 

For completeness we give also the following theorem on the existence of non- 
asymptotically stable T-periodic solutions for (|l.ip . In the theorem below, d{F,V) 
denotes the Brouwer topological degree of the vector field F e C"(R'^,M'^') on the 
open and bounded set ^ C M'' (see [H Ch. 2, § 5.2]). 

Theorem 2.9. Let g e C"'(]RxR'= x [0,1], R*^). A ssume that there exists an open 
hounded set V C M'"' such that go{v) 7^ for any v G dV and 
(vii) d{-go,V) < 0. 

Then there exists eo > such that for any e G (0, £0] system il.l]) has at least one 
non-asymptotically stable T -periodic solutions with Xe(0) G V. 

Proof. Since gQ{v) ^ for any v G dV then from Mawhin's Theorem (or [351 
Section 5]) we have that there exists £q > Q such that 

(2.5) d{-go.V)^d{I^x{T,-,e),V) for any e G (0, eo]- 

By p51 Theorem 9.6] for any asymptotically stable T-periodic solution of (|l.ip 
we have that d{I — x{T, ■,e),Bs{xe{0))) = 1 for (5 > sufficiently small. Therefore if 
all the possible T-periodic solutions of (jl.ip with e G (0, Eq] had been asymptotically 
stable, then the degree d{I — x{T, •, e), V) would have been nonnegative, contradicting 
(|vii)) and (1131). □ 

Remark 2.10. Assumptions Im]] and imply that d{—go,V) = 1 (see ]23\ 
Theorem 5.16]). 

Finally thinking in the application to the nonsmooth van der Pol oscillator, we 
formulate the following theorem which combines Mawhin's Theorem (see [31] or [32l 
Theorem 3], Theorem 12.51 and Theorem 12.91 In this theorem ([5o]i)(j) stays for the 
derivative of the i-th component of the function go with respect to the j-th variable. 

Theorem 2.11. Let g e C°(Rx il x [0, 1],M^). Let vq e be such a point that 
goivo) — and go is continuously differentiable in a neighborhood of vo. 

(a) If det (go)' (vq) 7^ then there exists Eq > such that for any e G (0,eo] 
system il.l]) has at least one T-periodic solution x^ such that x^iO) — > vq as 
£ ^ 0. 

(b) // 0j and (v\) hold and 

(2.6) det(go)'(wo) >0 and ([.9o]i)'(i)(wo) + ([.9o]2)'(2)(«o) < 0, 

then there exists Eq > such that for any e G (0,£o] system 11. V) has exactly 
one T-periodic solution x^ such that x^i^) — > wo as £ ^ 0. Moreover the 
solution Xe is asymptotically stable. 

(c) If det (go)' (vq) < 0, then there exists £0 > such that for any e G (0,£o] 
system il.l]) has at least one non-asymptotically stable T-periodic solution 
Xg such that Xe{0) — > wo as £ — ^ 0. 

Proof. Statement (a) is added for the completeness of the formulation of Theo- 
rem dTTTl and it follows from Mawhin's Theorem (see [3T| or [351 Theorem 3]). 

On the other hand it is a simple calculation to show that (|2.6p implies that all 
the eigenvalues of (goYivo) have negative real part. Therefore, assumption (jvl| of 
Theorem 12.51 is also satisfied and statement (b) follows from this theorem. 

Statement (c) follows from Theorem 12.91 Indeed since dct{goy{vo) < implies 
(see [331 Theorem 5.9]) that d{go, Bp{vo)) is defined for any p > sufficiently small 
and that d{go,Bp{vo)) = det(5o)'(wo) < 0. □ 
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3. Application to the nonsmooth van der Pol oscillator. In his paper [TH] 
Hogan first demonstrated the existence of a limit cycle for the nonsmooth van der Pol 
equation ii + e{\u\ — l)u + m = which governs the circuit drawn at Fig. 11.11 with the 
lamp characteristic i{u) = 5*0 + 5*1^ — 5'2w|w| whose derivative i'{u) = Si ~ 25*2 1?^! is 
nondifferentiable (see Nayfeh-Mook [35], §3.3.4, where the same stiffness characteristic 
appears in mechanics). In this paper we extend this study considering the van der 
Pol problem on the location of stable and unstable periodic solutions of the perturbed 
equation 

(3.1) u + e {\u\ — 1) it + (1 + ae)u = eXsint, 

where a is a detuning parameter and eAsini is an external force. We discuss with 
respect to the parameters a and A, under the assumption that e > is sufficiently 
small. 

Levinson's change of variables (see [28j, passing from Eq 2.0 to Eq 2.1) allows 
to rewrite equation (j3.1|) in a smooth form where the second Bogolubov's theorem is 
applied. But we remind that the idea of considering this example is to see what is the 
issues of the direct applying of Theorem 12.51 in comparison with the smooth results 
by Andronov and Vitt. 

Some function w is a solution of (|3.ip if and only if (zi, Z2) = {u, ii) is a solution 
of the system 

(3 2) ^ 

i2 = — zi + e[— azi — (|zi| — l)z2 + Asini]. 

After the change of variables 

\ _ f sill* \ f 2;i(t) 

Z2{t) I \ — sini cost J \ X2{t) 



system ()3.2p takes the form 
(3.3) 



ii = esm{—t)[—a{xiC0st + X2Smt) — 

— {\xi cost + X2 smt\ — 1) {—xi sint + X2 cost) + Asini] , 
X2 — ecos(— <) [—a(a::i cost + a:;2 sin i)— 

— {\xi cost + xi sini| — 1) {~X\ sint + X2 cost) + A sint] . 



The corresponding averaging function 50 7 calculated according to the formula (|1.2p . 
is given by 

^o.N [<7o]i(M,iV) = TTflTV - ttA + ttM - |MVM2TiV2, 

^ ' {gQ\^{M,N) = -vraA/ + TriV - |iVVM2 + 7V2, 

and it is continuously differentiable in M2\^|o}.. 

In short, by statement (a) of Theorem l2.11[ the zeros (A/, iV) G of this function 
with the property that det (go)'(Af, -^) 7^ 0, determine the 27r-periodic solutions of 
p.2p emanating from the solution of the unperturbed system 



, , wi(t) = Af cost + sint, 

^ ' U2(t) = -M sint + N cost. 

One has the following expression for the determinant 

(3.6) det(go)'(M,iV) = tt^{1 + a^) + ^{M^ + N^) - iny/ApTN^. 
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Following Andronov and Witt [5] we are concerned with the dependence of the am- 
plitude of the solution (|3.5p with respect to a and A, thus we decompose this solution 
as follows 

(3.7) ui{t) = Asui{t + (f)), U2{t) = Acos{t + (t)), 
where (M, N) is related to {A, (f) by 

(3.8) M = A sin (/), N ^ A cos (f). 
Substituting p.Sp into p.4p and (|3.6|) we obtain 

, . [go((^sin0,ylcos(/)))]i = —(4/3) • A\A\ sin0 + ttoAcos^ + TrAsin^ — ttA, 
^ ■ [5o((^sin0, Acos(/)))]2 = -(4/3) • cos(/)- 7raAsin0 + 7rAcos(/) 

and, respectively 

32 

(3.10) det (5o)'((^sin0, Acos(/.)) = 71^(1 + 0^) + - 2t:\A\. 

Looking for the zeros [A, (j)) of (|3.9p . we find the following implicit formula for deter- 
mining A: 




(3.11) A' \^a' 

Observe that the number of positive zeros of equation (|3.1ip coincides with the number 
of zeros of the equation A^ (^a^ + (l — -J^A)^^ = A^. To estimate this number we 
define 

and we have 

f'iA) - 





Since /' has one or two zeros then equation p. lip has one, two or three positive 
solutions A for any fixed a and A. In order to understand the different situations 
that can appear, we follow Andronov and Witt who suggested in [5] (see also |3j) 
to construct the so called resonance curves, namely the curves of dependence of A 
on a, for fixed A. Formula (|3.1ip is the equation of this curve. Some curves (|3.1ip 
corresponding to different values of A are drawn in Figure [SH] The way for describing 
these resonance curves p.lip is borrowed from [29l Ch. 1, § 5], where the classical 
van der Pol equation is considered. 

When A = the curve (|3.1ip is formed by the axis A — and the isolated 
point (0, 37r/4). When A > but sufhciently small the resonance curve consists of two 
branches: instead of A = we have the curve of the type / — / and instead of the point 
(0, 37r/4) we obtain an oval /' — /' surrounding this point. When A > increases, the 
oval /' — /' and the branch / — / tend to each other and, for a certain A there exists 
only one branch // — // with a double point P. The value of this A can be obtained 
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assuming that equation (|3.1ip has for a = a double root and, therefore, p.lOp 
should be zero. Solving jointly (|3.1ip and (|3.10p with a = we obtain A — 37r/16 
and P = 27r/8. If A > 37r/16 then we have curves of the type /// which take form V 
when A > crosses the value A = 9-\/37r/64. From here, if A < 37r/16, then equation 
(|3.1ip has three real roots when \a\ is sufficiently small, and only one root when \a\ 
is greater than a certain number which depends on A. When 37r/16 < A < 9A/37r/64 
we will have one, three or one solution according to whether a < ai, ai < a < 02 or 
a > a2, where 01,02 depend on A. The amplitude curves of type V provide exactly 
one solution of (|3.1ip for any value of a. The value A — 9'y37i'/64, that separates the 
curves where p. lip has three solutions from the curves where p. lip has one solution, 
is obtained from the property that (|3.1ip with this A has a double root for some a and 
thus this value of a vanishes (|3.10p . Therefore A — 9\/37r/64 is the point separating 
the interval (0, A) where the system formed by p. lip and 

39 

(3.12) ■k^(\ + c?) + —A^-2-k\A\=Q 

9 

has at least one solution from the interval (A, cxd) where p.lip - p.l2p has no solutions. 

In short we have studied the amplitudes of the 27r-periodic solutions of system 
(13. 2p depending on a and A. Whether a physical system described by (13. 2p possesses 
27r-periodic oscillations corresponding to some of these 27r-periodic solutions depend 
on whether some of these 27r-periodic solutions are asymptotically stable. To find the 
answer we use statement (b) of Theorem 12.111 Assumption (ji]) is obviously satisfied 
with = . Next statement shows that the right hand side of system p.3p satisfies 
©■ 

Proposition 3.1. Let vq e W- , vq ^ 0. Then the right hand side of itff.^) 
satisfies ^ for any a, A e M. 

The proof of the proposition is given in an appendix after this section. 

Thus we have to study the signs of p.lOp and ([go]i)^/(^sin0, ^cosc/)) 
+ {[g o]2y jsf { A sin (f), A cos (j)). We have 

(Mi)m(M,7V) + {[goh)'NiM,N) = 2 (n - 2^ AP + N^) , 

and therefore the conditions for the asymptotic stability of the 27r-periodic solutions 
of ([321) near are 



(3.13) TT'^{l + a^) + —{M"^ + N^) -iny^Ap'+Jp > 0, 
and 

(3.14) 2(n~2VM^ + N^^ < 0. 

Substituting (13. 8p into the inequalities p.l3p and p.l4p . we obtain the following 
equivalent inequalities in terms of the amplitude A 

(3.15) TT^{l + a^) + ^A^ -27r\A\>0 

9 

and 



(3.16) 



27r-4|A| < 0. 
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Conditions (|3.15p and (|3.16p mean that the asymptotically stable 27r-periodic solu- 
tions of (|3.2p correspond to those parts of resonance curves under consideration which 
are outside of the ellipse (|3.12p and above the line A = tt/2. All the results are col- 
lected in Figure 13.11 from where it is easy to see that for any detuning parameter a 
and any amplitude A > 0, equation p.ip possesses at least one asymptotically stable 
27r-periodic solution with amplitude close to A obtained from (|3.1ip . Among all the 
asymptotically stable 27r-periodic solutions of p.ip . there exists exactly one whose 
fixed neighborhood does not contain any non-asymptotically stable 27r-periodic solu- 
tion of (|3.ip for sufficiently small e > 0. The amplitude of this asymptotically stable 
27r-periodic solution is obtained from (|3.15p - (|3.16p . 



■A 









J. 3 
















jr 
















/ U l-E 














1 » 


= '' I'i ' 


















1 






















¥.---''' 




' ' 










_ - 












■7-' 














"T 


-1,5 


-a 










1 


IS 



Fig. 3.1. Dependence of the amplitude of stable (solid curves) and unstable (dash curves) 2tt- 
periodic solutions of the nonsmooth periodically perturbed van der Pol equation i3.Hl on the detuning 
parameter a obtained over formulas i3.lH) . {3715]) and 13. 16V for different values of X. The curve I 
is plotted with A = 0.4, II with X = 37r/16, /// with some X = ^^OA e (37r/16, 9V37r/64), IV with 
X = 9^371/64, V with X = 1.5. Point P is 2/^3. 



To compare the changes due to nonsmoothncss in the behavior of the resonance 
curves, we give in Figure lHT^ the resonance curves of the classical van der Pol oscillator 



(3.17) 



il + elu^ -1] 



(1 + ae)u = eX sint, 



which can be found in [2, Fig. 4] or in Ch. I, § 16, Fig. 15]. 

The key formulas for Figure 13.11 can be also easily comparable with the key formulas 
for Figure In fact the corresponding expressions (|3.1ip - p.l2p and (|3.13p - (|3.14p 
are (see the formulas (5.21)-(5.22) and (16.6)-(16.7) from j29j) 




A^ 



16 



-A^ - 0, 
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Fig. 3.2. Dependence of the amplitude of stable (solid curves) and unstable (dash curves) 2it- 
periodic solutions of the classical periodically perturbed van der Pol equation [S-IT^ on the detuning 
parameter a for different values of X. Following Andronov-Witt (see f^, Fig. 4) curve I is plotted 
with X = VOA, II with X = 4v^/9, /// with some 4v^/9 < A < ^32/27, IV with X = V'32/27, V 
with A = 2. Point P is 2/\/3. 



and 

l + a^~ (M^ + N^) + -^(M^ + N^f > 0, 
16 

2 - (M^ + N^) < 0, 

respectively, when considering the classical van der Pol equation p.l7p . 

It can be checked that the eigenvectors of the matrix (go)' {{A sin (p, A cos (j))) are 
orthogonal only for A = that is Theorem 3.6 from Buica-Daniilidis paper [S] can 
not be applied. At the same time, assumption (H2) from 9^ is not satisfied for our 
problem (see Remark 12.81) . 

4. Appendix. Proof of Provosition [X71 As before, [v\i is the i-th component 
of the vector v E M^. Let g{t, v) = |[t;]i cost + [v\2 s\nt\ and notice that it is enough 
to prove that g : [0, 27r] x ^ M satisfies (jvj). In the case that [wo]2 7^ 0, denote 
9{v) = arctan(— [w]i/[w]2), while when [wo]2 = 0, denote 6{v) = arctan(— [w]i/[?;]2) 
for [wo]ib]2 < 0, 9{vq) — 7r/2 and, respectively, 6{v) = arctan(— [u]i/[w]2) + tt for 
bo]i[w]2 > 0. In any case notice that the function v ^ 9{v) is continuous in every 
sufficiently small neighborhood of wq- Fix 7 > 0. Let M be the union of two intervals 
centered in 9{vq) (when 9{vq) < 0, take 0{vo) + 2n instead) and, respectively, 9{vo)+tt, 
each of length 7/2. Denote them Mi and M2. Take S > such that 9{v) G Mi for 
all V G Bg{vo). Of course, also 9{v) + tt G M2 for all |jt; — vo\\ < 5. This implies that 
for fixed t G [0, 2tt] \ M, [v]i cost + [v]2 sint has constant sign for all v e Bg{vQ), that, 
further, gives that g{t,-) is differentiable and g[,{t,v) = gy{t,vo) for all v £ B-^{vo). 
Hence (jvj) is fulfilled. □ 
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